Newtonian Perturbations on Models with Matter Creation 
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Creation of Cold Dark Matter (CCDM) can macroscopically be described by a negative pressure, 
and, therefore, the mechanism is capable to accelerate the Universe, without the need of an ad- 
ditional dark energy component. In this framework we discuss the evolution of perturbations by 
considering a Neo-Newtonian approach where, unlike in the standard Newtonian cosmology, the 
fluid pressure is taken into account even in the homogeneous and isotropic background equations 
(Lima, Zanchin and Brandenberger, MNRAS 291, LI, 1997). The evolution of the density contrast 
is calculated in the linear approximation and compared to the one predicted by the ACDM model. 
The difference between the CCDM and ACDM predictions at the perturbative level is quantified by 
using three different statistical methods, namely: a simple x 2 - anai y s i s m the relevant space param- 
eter, a Bayesian statistical inference, and, finally, a Kolmogorov-Smirnov test. We find that under 
certain circumstances the CCDM scenario analyzed here predicts an overall dynamics (including 
Hubble flow and matter fluctuation field) which fully recovers that of the traditional cosmic con- 
cordance model. Our basic conclusion is that such a reduction of the dark sector provides a viable 
alternative description to the accelerating ACDM cosmology. 



I. INTRODUCTION 

The growing evidences for cosmic acceleration as, e.g., 
the Supernova la observations [lHj], is challenging cos- 
mologists and theoretical physicists. The neediness of 
some new ingredient in the cosmic recipe, in order to pre- 
serve Einstein's Equations, inspired both communities to 
conservatively invoking the simplest available hypothesis, 
namely, a cosmological constant, A. 

Nevertheless, the identification of A with the quan- 
tum vacuum has brought another problem: the estimate 
of theoretical physicists that the vacuum energy density 
should be 120 orders of magnitude bigger than the mea- 
sured A value. This is the "old" cosmological constant 
problem (CCP) gj]. The "new" problem [6| asks why is 
the vacuum density so similar to the matter density just 
now? Many solutions to both problems have been pro- 
posed in the literature @, Q ■ The majority of them are 
plagued with no physical basis and/or many parameters. 

On the other hand, some authors have also investi- 
gated a class of models where the creation of cold dark 
matter may result on a pressure which is negative (named 
CCDM, Creation of Cold Dark Matter), thereby provid- 
ing a mechanism for cosmic acceleration. In the current 
literature, the leitmotiv of such models is to reduce the 
dark sector (dark energy + dark matter) in the frame- 
work of general relativity since the existence of dark en- 
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ergy seems to be less secure than dark matter [§, [T(| • The 
extra bonus is to solve the coincidence and cosmological 
constant problems. 

The search for dark matter accelerating models started 
almost one decade before the SNe la observations. Ini- 
tially, Prigogine and coworkers [ll| argued that the 
gravitationally-induced particle creation could consis- 
tently be discussed in the realm of the relativistic non- 
equilibrium thermodynamics. Later on, their macro- 
scopic formulation was clarified by Calvao, Lima and 
Waga through a manifestly covariant formulation pjj]. 
The inclusion of the back reaction in the Einstein Field 
Equations (EFE) via an effective pressure (which is neg- 
ative for an expanding space-time) opened the way for 
cosmological applications. As a consequence, several in- 
teresting features of cosmologies where the dark sector is 
reduced due to the creation of CDM matter have been 
discussed in the last decade [HI, H3] . 

More recently, a new CCDM cosmology was proposed 
by Lima, Jesus and Oliveira (from now on LJO model) 
which mimics the Hubble expansion history of the ACDM 
model, at least at the level of the background equations 
.15]. 

The quoted LJO model is quite interesting, once that it 
perfectly mimics the ACDM cosmic history, thereby also 
providing a good fit to current cosmological data (SNIa, 
BAO and CMB shift parameter, total age, etc) with the 
same number of parameters of the ACDM, but without 
the CCP. In this sort of models, the quantum vacuum 
energy is canceled out by some not yet known physical 
mechanism. A basic advantage of this kind of scenario 
comes from the fact that a vacuum energy density which 
is null or negligible is more acceptable than a vacuum 
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energy density finely tuned like in the ACDM model. 

Naturally, although providing good fits to background 
data, it would be interesting to investigate to which 
level this similarity is preserved, for instance, by ana- 
lyzing the evolution of small fluctuations predicted by 
the CCDM equations. Preliminary studies on this sub- 
ject has already been accomplished by Basilakos and 
Lima [16| . however, focusing on the theoretical conse- 
quences to cluster abundances at different redshifts (see 
also [TtJ )• By using the Press-Schechter formalism, they 
also investigated the cluster-size halo redshift distribu- 
tion by confronting the results with future cluster sur- 
veys (eROSITA satellite and Sunayev-Zeldovich survey 
based on the South Pole Telescope). We would like to 
stress that in the latter papers we ignored possible con- 
tributions from pressure at the perturbative level. 

In the present paper, we are basically interested to 
determine whether the LJO model provides a realistic 
description at the perturbative level, however, by fit- 
ting the available fluctuation data, like the growth rate 
of clustering. The difference between the CCDM and 
ACDM predictions at the perturbative level is quantified 
using tree different statistical methods, namely: a simple 
X 2 -analysis in the relevant space parameter, a Bayesian 
statistical inference, and, finally, a Kolmogorov-Smirnov 
test. As we shall see, the CCDM reduction dark sector 
provides a viable alternative description to the accelerat- 
ing ACDM cosmology. 

The work is planned as follows. The general relativis- 
tic approach for CCDM models is presented in section 
II. In sections III we show how the Neo-Newtonian treat- 
ment adopted here recovers the relativistic formulation 
whereas, in section IV, we derive the associated evolu- 
tion equation for the density contrast. In section V, we 
discuss in detail the linear growth factor of matter per- 
turbations, the CCDM theoretical predictions regarding 
the evolution of the growth rate of clustering, and the 
corresponding statistical analyzes. The main conclusions 
are summarized in VI. 
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where T is the creation rate of CDM particles and H = 
a I a is the Hubble parameter. 

It is readily checked that equations ([TJ and @ lead to 
the continuity and acceleration equations: 
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The above Eqs. dU[S|) imply that the model is fully 
determined by the creation rate T, or more precisely, by 
the ratio T/H. If V << H, the particle creation pro- 
cess is negligible, and the standard dust filled models are 
recovered. In the LJO scenario, the phenomenological 
creation rate of cold dark matter has been defined by: 
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where Cl\ (termed a in the original LJO notation (l5j ) 
is a constant parameter, p C Q is the present day value of 
the critical density, and the factor 3 has been added for 
mathematical convenience. 

Now, by inserting the above expression into the energy 
conservation as given by Q one obtains 



p + 3Hp = Tp = 3fl A p c0 H, 



(7) 



which can be readily integrated to give a solution for the 
energy density 



P=(pQ- &APco)a 3 + VLaPc 



(8) 



II. CCDM COSMOLOGY IN THE MATTER 
ERA: RELATIVISTIC FORMULATION 

The background cosmological equations of the model 
with creation of cold dark matter (CCDM) have the fol- 
lowing form (for simplicity, we are neglecting the contri- 
butions of the radiation and baryonic components): 
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where p is the CDM density, p c is the creation pressure, 
a(t) is the scale factor, and an overdot means time deriva- 
tive. In the case of constant specific entropy (per parti- 
cle) .the creation pressure for a CDM component is given 
by PtH 



In terms of the cosmic history, the equivalence between 
this model and the ACDM model can be seen directly 
through the evolution equation of the scale factor func- 
tion. As one may check, by inserting the expression of 
the creation pressure p c in Eq. @ we obtain: 
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which should be compared to: 



2aa + a 2 + K - Aa 2 = 0, 



(9) 



(10) 



provided by the ACDM model. The above equations im- 
ply that the models will have the same dynamical be- 
havior when we identify the creation parameter by the 
expression = A/3i?o 2 = ^A- Further, the factor 
f2 m = Vl m — f^A) where Vl m = can also be identified 

PcO 

as an 'effective' CDM density parameter, in comparison 
to the ACDM model EE!. 
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In the context of a spatially flat geometry (K = and 
&m = l)i one may show that the inflection point, that is, 
the point where the Hubble expansion changes from the 
decelerating to the accelerating regime (a(tj) — 0) takes 
place at: 

r - -1 1/3 

a/ = . (11) 

|_2fi A J 

As an example for fi m = 1 — Cl\ = 0.30 we find a/ ~ 0.60 
(or zj ~ 0.67). Finally, the necessary condition for an 
inflection point in our past is ai < 1, which leads to the 
condition -| < ft a < 1- 



linear perturbation is equivalent to the full general rel- 
ativistic formulation, at least when some conditions are 
imposed [2lJ . Indeed, applications of this Neo-Newtonian 
approximation are not restricted to nonrelativistic mat- 
ter, and the high accuracy of the approximation has also 
been proved for different epochs and even for scales larger 
than Hubble radius, as well as, for the spherical gravita- 
tional collapse (22), [23[. Therefore, it is reasonable to 
expect that the Neo-Newtonian approach will provide a 
good approximation for the evolution of small density 
fluctuations in the CCDM case. 



A. Recovering the Relativistic CCDM Model from 
the Neo-Newtonian Approach 



III. NEO-NEWTONIAN FORMULATION FOR 
CCDM COSMOLOGIES 

The relativistic dynamics of CCDM cosmologies is the 
same of a single fluid with density p and pressure p c (see 
Eqs. ([T][5])). In addition, whether the equivalent fluid has 
an equation of state (hereafter EoS) parameter, w = p/p, 
taking a look at Eq. ([3]), we may see that this fluid 
(p = Pc), in general, has a time varying equation of state 
parameter, w = —jjj- 

In this connection, it is interesting to notice that flu- 
ids with non-vanishing pressure were also studied by 
Lima, Zanchin and Brandenberger [l8j within a 'quasi'- 
Newtonian perspective. Follo wing earlier ideas developed 
by McCrea [19| and Harrison [20|, they proposed the so 
called Neo-Newtonian description in order to reproduce 
the relativistic equations with pressure either at the back- 
ground and perturbative levels. The basic equations of 
such a description are given by: 

(£) + (u.V r )u = -V r *-^ (12) 

V dt J r P+Pc 

(It) + V r- {pu)+ Pc \7 r -u = (13) 

= 4ttG(p + 3 Pc ) (14) 

where we have already included the creation pressure, 
p c . These equations are named Euler, continuity and 
Poisson equations, respectively. The modified continuity 
equation is due to Lima et al. [HI, in order to account 
correctly to pressure effects. As one may check, such 
equations reproduce the FRW type equations with pres- 
sure in the homogeneous and isotropic case (p c = p c (t), 
p = pit)), and can also be consistently perturbed for any 
given equation of state (see discussion below) . In partic- 
ular, they also shown that the linear perturbed version 
of these equations yield the correct growing mode for the 
density contrast of a fluid with p = wp, with constant 
w. For the case that w is a time dependent quantity, the 



Let us now discuss how the basic relativistic equations 
of CCDM model can be recovered by the set of equa- 
tions (TT2"]) , (fT"3]) . and (|14p . In order to show that, we 
first remark that homogeneity and spatial isotropy of the 
unperturbed model implies u = Hr, according to the 
Hubble's law. Then, Eq. (flU)) can be rewritten as: 

p + 3Hp = Tp. (15) 

Assuming the creation rate T as in the Eq. © , the above 
equation can readily be integrated, to give exactly the 
same expression for the density evolution in Eq. (jS| . The 
Euler equation (TT2"]) can be rewritten as 

Hr + Hr—{Hr) = (H + H 2 )r = -V r <5>. (16) 
or \ J 

While the integration of the Poisson equation (p~4|) yields: 

V r $ = 47rG(p + 3p c )^, (17) 

Finally, by combining Eqs. (1161) and (1171) . we derive 
the acceleration equation: 

a 4nG , „ . 

- = ^(P + 3 PC ), (18) 

where, as in the relativistic case, there is the presence of 
the creation pressure term, which can give rise for accel- 
eration since it is negative for the expanding Universe. 
Using the expressions (|3]) and ©, the equation above 
can be integrated resulting in the following expression: 

"\K 

8ttG P = 3H 2 + — . (19) 

a z 

where K is an arbitrary integration constant. This is 
Friedmann equation (TTJ). 

Thus, Eqs. (fT8|) and (fT9| show that the background 
cosmological equations can effectively be recovered by the 
Neo-Newtonian formulation. This result is valid regard- 
less of the specific form assumed to the creation rate T, 
and should be compared with the incomplete treatment 
adopted by Roany and Pacheco [Hj], and, critically, re- 
discussed by Lima et al. [Hj]. 



4 



IV. NEO-NEWTONIAN DENSITY 
PERTURBATIONS 

In general, a perturbative analysis in cosmology re- 
quires a full relativistic description since the standard 
non-relativistic (Newtonian) approach works well only 
when the scale of perturbation is much less than the Hub- 
ble radius and the velocity of peculiar motions are small 
in comparison to the Hubble flow (2(|. However, as re- 
marked earlier, such difficulties have been circumvented 
by the Neo-Newtonian approximation adopted here, and, 
therefore, such an approach may provide a good feeling 
about the behavior of the density perturbations. 

To begin with, let us transform the Neo-Newtonian 
Equations (TT^|) - p^)) to comoving coordinates x, which 
are related to the proper coordinates r by jig, [271] : 



r 

x = — . 

a 



(20) 



In these coordinates, the basic quantities u, p and $ can 
be rewritten as: 



u = ax 



ax = ax 



v. 



(21) 



' 1 *i 

p—+3H(6 Pc 



PcS) 



P + Pc 



^v-V<5 = 0, (28) 



V 2 </> = 4irGa 2 {p5 + 35p c ). 



(29) 



It should be stressed that the above Eqs. (|2"T)) - (f2"9"|) are 
obtained only by using the unperturbed background Eqs. 
(El)- (El) - In order to simplify the above equations it is 
suitable to dehne the quantities: 
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Sp c 
5p 



Sp c 
pS 



(30) 



Now, Eqs. (|2"71) - (|2"9")) can be rewritten as: 
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~dt 
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Hv 



(v • V)v 



(l + <* 



eff 
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(31) 



+3H(cl fr 



w)S- 



l + w + (l + c 2 eff )S 



V-v- 



(32) 



p = p(t)[l + 5(x,t)] 



(22) 



V 2 ^ = 4 7 rGa 2 p ( 5(l + 3cL f ). 



(33) 



Pc =Pc(t) +5p c (x,t), 



2-kGo 2 o _ 2 
$ = <f>+——(p + 3p c )x 2 , 



(23) 



(24) 



where v is the velocity perturbation (peculiar velocity), 
6 is the dark matter density contrast, dp c is the pertur- 
bation of creation pressure and </> is the peculiar grav- 
itational potential, which generates the peculiar accel- 
eration. Here, variables with an overbar represent the 
background quantities. 

Following standard lines, we change the proper coordi- 
nates (r, t) to comoving coordinates (x, t) and using the 
differential operator identities 



V x — V = aV r 



(25) 



These equations can already be compared to the equa- 
tions of [23j], where they study the perturbations of a 
mixture of non-interacting fluids using various EoS pa- 
rameters. These are identical to their equations for the 
case of one fluid with an EoS parameter w. 

Here we are interested only on the linear order of the 
perturbations. In this case, we see that Eq. (l33t is al- 
ready linear while Eqs. (|3"Tj) - (p?2"|) are reduced to: 
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V0 
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-V • v = 0. 



(34) 



(35) 



and 



Now, by calculating the divergent of Eq. (IM)) and insert- 
ing (I33p in the resulting equation, it is easy to see that 
Eq. d33} takes the form 
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dt 
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Eqs. (rr2^ - flT4"|) can be rewritten as: 



dv a (v ■ V)v VA 

— + -v + ^ >— = — r 

dt a a a 



V5pc 



a{p + Pc)' 



(27) 



where we have defined 9 = V • v and assumed 



'eff 



"eff 



(t) and that the spatial dependence of 6 is propor- 
tional to e . Now, we may isolate 9 = V • v on (|35[) . 
and replace it into (|3"61 to finally find: 
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dt 2 



H{2 + 3c 2 eff - 3w) 



3{H + 2H 2 )(c 2 eff - w) + 3H 



C e//-( 1 + C e//) T 



■w 



4irGp(l + w)(l + 3c 2 
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1.2 2 
K °eff 



86 
dt 



eff' 



6 = 0, (37) 



which is the same result found by 21 1, in the context of case, w = —jfi, we find: 
a fluid with EoS parameter w. Now, recalling that in our 
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where we have used Eq. (|5|). For numerical purposes it variable, r\ = ln(a(f)): 
is interesting to write (J39J) in terms of a new independent 
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where the prime denotes derivative with respect to 77. If and (|40p simplifies to: 
we live in a spatially flat Universe, K = 0, then = |, 
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Now, if the pressure perturbation vanishes 16j, c; 
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V. CASE STUDY: LJO MODEL 

The creation rate of the LJO model reads (see Eq. ©) 



r = 30 A [ ^ ) //. 



whereas the Hubble parameter is given by: 



H = Hn 



(O m - tt A )a 3 + tt A + (1 - Q m )a~ 



(43) 
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Po 



In the flat case, O m = 1, and the 



Hubble parameter reduces to: 



H = H 



(i - n A ) a - 3 + fi/ 



1/2 



(45) 



For the above creation rate, the linear order perturba- 
tion equation (|40l) can be rewritten as: 



5" + F(rj)8' + G(rj)S = 0, 



(46) 



' where the functions F(rj) and G(rj) are given, in the flat 

(44) case (g5]), by: 
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f4S) 



In what follows, the evolution of the density contrast 
as given above will be numerically solved by assuming the 
same initial conditions of the Einstein - de Sitter grow- 
ing model, namely: <5(a,;) = a, and <5'(a,;) = 1, where 
dj = 1CP 3 . It is natural to impose such conditions be- 
cause the particle production in the LJO model is neg- 
ligible at high redshifts with the model reducing to the 
standard dust filled FRW cosmology. In addition, we also 
remark that the integration requires a choice of c 2 ^ (ap- 
pearing in the functions F(rf) and G[rf)) or, equivalently, 
the perturbation of the creation pressure. 

To begin with, let us recall that Sp c is a new degree 
of freedom [2^, [28| which does not depend explicitly on 
the known fluid variables, like S. Once the pressure p c 
form is given, one may try to estimate Sp c . For the LJO 
model, by using Eqs. §5$ and (|4"3")l one finds: 

p c = -Q A p c0 , (49) 

which is time independent. Thus, one could expect Sp c = 
0, which would correspond to c 2 ^j = 0. Inserting this 
condition into (|4"?| and (|48p we obtain: 

(1 - fi A ) + l&n A e 3 ^ 



2(1 - A + fl A e 3r >) 



(50) 



G(v) = 



9(1 -0 A ) 2 



3(5f2 A e 3? ? -2 + 20 



2(1 - o A + n A e 3r >) 2 i-n A + n A e 3r i 



This case is also interesting because the condition for 
the Neo-Newtonian perturbed equations to be equivalent 
to the full general-relativistic treatment for a single fluid 



is that c 2 



eff 



o mi. 



Another interesting choice to the effective sound speed 
is Cgjj = c 2 = Pc/ p. In this case, the non-adiabatic pres- 
sure, Sp na d ~ Spc — c 2 Sp vanishes, and one has only adi- 
abatic perturbations as suggested by CMB observations. 
However, the sound speed of matter in the presence of 
particle production reads 



HY' - YH' - Y(3H - Y) 
3H(3H - Y) 1 



(52) 



which vanishes identically when the creation rate is given 
by Eq. (|4"3"]) since p c = 0. Therefore, in the LJO frame- 
work, this choice for the effective sound speed reduces 
to the case earlier analyzed (c 2 ^ =0). Of course, since 
c 2 ff is a new degree of freedom, a possibility is to choose 
it as a general time dependent variable, say, c 2 jj = w. 
The treatment involving c 2 ^ as a new degree of freedom 
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FIG. 1: 8 as a function of a for the case cf// =0. As discussed 
in the text, for LJO models this is equivalent to the case 
c^ff = <? a because the sound speed is zero. For comparison 
we also show the fiat ACDM case, with Q m = 1 — S1a = 0.3 
and Cgjj =0. 



is relevant when the unperturbed fluid equations evolve 



out of thermodynamic equilibrium as in the case of LJO 
model. 

In Figure 1 , we display the evolution of the density con- 
trast for LJO model as predicted by the Neo-Newtonian 
approach. Note that the increasing mode perturbation 
grows until its maximum value after which it decays in 
the course of the evolution. In particular, for il A = 0.2, 
we see that the predicted density contrast is indistin- 
guishable from that of ACDM cosmology. However, we 
stress that such a value is not favored by the background 
tests which prefer & A = 0.7 [151. [la]. 

In Figure 2, we display the case c 2 ^ — w — —F/3H. 
As shown there, the evolution of the density contrast has 
been obtained for different values of fl A . Evidently, for 
0.2 < Cl A < 0.4 the evolution is similar to the ACDM 
behavior. 

As it appears, the simplest possibility to the quantity 
c^ff is to consider it as a constant free parameter and 
find out which value is preferred from observations. In 
this case, c 2 ^' = and now we have: 



F(r>) 
G{v) 



(1 - 00(1 + 6c 2 eff ) + 2ft A e 3 "(8 + 



2(1- 
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H 2 {l~n A + n A e 3 *>) 



(53) 
(54) 
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FIG. 2: 8 as a function of a for the case c 



eff 



Also 



shown, for comparison, is the flat ACDM case, with f2 m = 
1 — Q.a = 0.3 and c^tt = 0. In this case we have considered 
the limit of large scales (k = 0). 

In Figure 3, we show the evolution of the density con- 
trast when Cgjj is taken to be a constant free parameter. 
The selected values of c 2 e j^ are also indicated. Again, we 
see that there is a range of values for which the CCDM 
matter fluctuation field mimics the behavior of the con- 
cordance ACDM model. 



A. The growth rate of clustering 

We would like to end this section with a discussion 
on the evolution of the well known indicator of cluster- 
ing, namely the growth rate 27J. This is an efficient 
parametrization of the linear matter fluctuations 5(a) 
which has the following functional form: 



/(*) 



din6 
cflna 



-(1 + *) 



diaS 
dz 



(55) 



In Table I, we show the existing growth rate data with 
the corresponding error bars, the associated redshifts and 
related references. 

Let us now attempt to place constraints on the relevant 
parameters by performing a standard x 2 minimization 
procedure between the observationally measured growth 
rate and that predicted by the LJO cosmology. The best 
fit to the set of independent parameters (£l A ,c 2 e ff) can 
be estimated by using a \ 2 statistics with 



x 2 (n A ,c 2 eff ) = 



7 

E 

i=l 



fobs( z i) /model (^i i ^A; ^ef f) 



(56) 
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a 

FIG. 3: 5 as a function of a for the case in which is a free 
parameter. Also shown, for comparison, is the flat ACDM 
case, with Q m = 0.3 and c^// = 0. In this case we have 
considered the limit of large scales (k = 0), and CIa = 0.7. 

TABLE I: Data of the growth rate of clustering. The cor- 
respondence of the columns is as follows: redshift, observed 
growth rate and references. 



z 


fobs 


Refs. 


0.15 


0.49 ±0.10 


[29, 30] 


0.35 


0.70 ±0.18 


[31j 


0.55 


0.75 ±0.18 


[32] 


0.77 


0.91 ±0.36 


[30] 


1.40 


0.90 ±0.24 


[33] 


2.42 


0.74 ±0.24 


m 


3.00 


1.46 ±0.29 


[35] 




z 



FIG. 4: Upper Panel: Likelihood contours in the (Qa,^//) 
plane. The contours are plotted where — 21n(£/£ max ) is equal 
to 2.32, 6.16 and 11.83, corresponding to la, 2a and 3a con- 
fidence level. The perpendicular dashed line separates the 
regions in which the flat LJO model can accommodate an ac- 
celerated expansion (Qa > 1/3), equivalent to the standard 
ACDM cosmology. Bottom Panel: Comparison of the ob- 
served (solid circles see Table I) and theoretical evolution of 
the growth rate of clustering f(z). Note that we have imposed 
Qa = 0.7 as suggested by the background tests [TBI. [T^]. The 
lines correspond to the LJO (solid curve with C^ff = —1) and 
the ACDM (dot-dashed curve with Qa = &a and c^/< = 0) 
models. 



where ct; is the observed growth rate uncertainty. Note 
that we sample the free parameters as follows: Q\ £ 
[0.1, 1] and c 2 eff £ [-1.3, 2] in steps of 0.001. 

In Figure 4 (upper panel), we show the la, 2a and 3<r 
confidence levels in the (Cl\,c^ff) plane. It is evident 
that Cgjj is degenerate with respect to (Ia and that all 
the values on the interval — 1.3 < c 2 e ^ < 2 are acceptable 
within the la uncertainty. As one may check, the c 2 ^ 

parameter as a function of Oa is well fitted by a power 
law having the form: 

cl ff = 0.333(±0.011)^°- 989(±0 ° 21) - 1-5 ■ 

It is also evident that the likelihood analysis puts some 
constraints on the value of c 2 ^ once including the nec- 
essary condition for an inflection point, that is, a tran- 
sition from a decelerating to an accelerating regime in 
our past. Thus, for (Ia > 1/3, (see dashed line in the 
upper panel of figure 4) we find that c 2 ^ lies in the in- 
terval — 1.24 < c 2 ^ < —0.72. In this framework, if wc 
marginalize over J7a = 0.7, then the likelihood analysis 
provides a best fit value of c 2 ,^ ~ — 1 with \ 2 — 4.03 
(dof = 6). The latter result remains unaltered within a 



physical range of 0.65 < fl\ < 0.74. Of course, for com- 
parison we perform the same analysis also for the ACDM 
model (fi A = 0.7) and we find X 2 ^ 4.88 (dof = 7). 

In the bottom panel of Figure 4, we present the LJO 
growth rate of clustering (solid line). We find that for 
z < 1 the LJO growth rate is somewhat less with respect 
to that of the usual A cosmology (dot-dashed line) and 
such a difference increases for extremely low redshifts. 
However, at relatively large redshifts (z > 1) the growth 
rate of the LJO model is indistinguishable from that of 
the ACDM model. This should be expected because the 
flat LJO model reduces to the Einstein - de Sitter cosmol- 
ogy at intermediate and high z's (negligible particle pro- 
duction). Note that the measured fobs(z) are presented 
in the bottom panel of Fig. 4 by the filled symbols. 

We further explore our statistical results by using a 
Bayesian statistics (see for example 36]), in which the 
corresponding estimator is defined as: BIC = x 2 +k n h\N 
[where k n is the number of free parameters and N{= 7) is 
the number of data points used in the fit] . The next step 
is to estimate the relative deviation between the two mod- 
els ABIC = BIC ljo — BIC a- In general a difference 
in BIC of ABIC > 6, is considered evidence against 
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that model which occurs the larger BIC. In our case, 
we find ABIC ~ 1.01 which implies that the LJO model 
with (^A,Cg^j) = (0.70, —1) fits very well the growth rate 
data. The latter result holds also for 0.65 < < 0.74. 

Secondly, we compare the growth rate of clustering be- 
tween data and models via a Kolmogorov-Smirnov (KS) 
statistical test respectively by computing the correspond- 
ing consistency between models and data (Vks)- Al- 
though both cosmological models fit well the data it 
seems that the LJO model provides a slightly better fit 
(Vks — 0.997) with respect to that of the usual A cos- 
mology (Vks — 0.883). Note that the KS test between 
the two cosmological models gives Vks — 1- 

Based on the above statistical tests it becomes evident 
that in the LJO model (0.65 < fl A < 0.74, c 2 eff = -1), 
the corresponding Hubble flow as well as the matter fluc- 
tuation field resembles that of the traditional ACDM cos- 
mology without the need of the required, in the classical 
cosmological models, dark energy. 

At this point one may ask about the possibility of fu- 
ture detection or at least what is a clear cosmic signature 
of this gravitationally induced creation process of cold 
dark matter particles. As it is widely known, some very 
massive dark matter particle candidates like the wimpzil- 
las (Mw ~ 10 13 GeV) can be copiously produced only 
at the very early stages of the Universe, mainly at the 
end of inflation [lOL [37J, and, as such, this kind of candi- 
date does not fit in our phenomenological treatment for 
continuous matter creation. In principle, a more realistic 
scenario is provided by the semi-classical approach pro- 
posed by Parker and collaborators [Hj], where massive 
particles (associated to a real scalar field) can be con- 
tinuously created during the expansion of the Universe. 
On the other hand, since the current dark matter detec- 
tion experiments rely on the physical properties of the 
dark matter (mass, cross section, etc), we are unable to 
identify (based only on our phenomenological approach) 
which is the specific candidate for the continuous cold 
dark matter production discussed here. However, by as- 
suming that the created mass are of the form of neu- 
tralinos (Mjv ~ lOOGeV), one may show that its present 
creation rate, T neu ~ 10~ n cm~ 3 yr -1 , has not changed 
appreciably in the last few billion years when the Uni- 
verse entered its accelerating phase (l5| . 

VI. CONCLUSIONS 

The problem related to the nature of non-baryonic 
component filling the observed Universe is usually re- 
ferred to as the non-baryonic dark matter problem. In 
the last three decades, many possible candidates from 
particle physics have been proposed to describe such a 
dark matter component [10]. In the framework of gen- 
eral relativity, it is also interesting to know whether the 
dark sector (dark matter + dark energy) could be prop- 
erly reduced only to the dark matter component with 
the conventional ACDM emerging as an effective descrip- 



tion [ll| [l6| . Naturally, unlike the decelerating dust-filled 
Einstein-de Sitter Universe, an extra mechanism must be 
responsible for the present observed accelerating stage. 

In this context, we have performed a Neo-Newtonian 
description of relativistic CCDM models, that is, mod- 
els endowed with gravitationally induced production of 
cold dark matter particles. The complete equivalence of 
the Neo-Newtonian cosmological approach with the gen- 
eral relativistic background equations regardless of the 
specific form adopted to the particle production rate, 
r, has been established. In the same vein, the general 
perturbed equations has also been derived in this frame- 
work. Due to the form of the negative creation pressure 
(p c = —Tp/3H) we have determined the general differ- 
ential equation to the contrast density by assuming a 
full equivalence to a fluid model with an EoS parameter 
w = —T/ZH. The resulting equations of section IV are 
also valid for arbitrary forms of T. 

In section V, we have focused our attention to a specific 
CCDM cosmology, namely, the LJO model |l5|. The 
interest for this special class of CCDM cosmology comes 
from the fact that its expanding history is equivalent to 
the standard ACDM model. 

Three different perturbed scenarios depending on the 
choice of the c 2 ^ were analyzed. In the first one c 2 e ^ — 
c 2 s = and the linear perturbation is equivalent to that 
predicted by ACDM only up to redshifts of the order 
of one. It was found that the increasing mode decays 
rapidly after z <~~> 1 with the redshift marking the begin- 
ning of attenuation depending on the parameter (see 
Fig. 1). We also show that the nonadiabatic case with 
c 2 jf — u is equivalent to a perturbed ACDM, however, 

only for 0.2 < < 0.4 (see Fig. 2). Finally, when c 2 ff 
was treated as a free constant parameter, the perturbed 
evolution is like ACDM for (l\ = 0.7 (see Fig. 3) which is 
exactly the same value preferred by the expanding cosmic 
history [15]. For all these models, it should be interest- 
ing to obtain the predicted matter power spectrum in 
the context of the relativistic formulation. This work is 
in progress and will be published elsewhere. 

For completeness, we also performed a detailed statisti- 
cal analysis based on the observed growth rate of cluster- 
ing in order to constrain the free parameters (Q,^,c 2 ^). 
Interestingly, when combined with the background tests 
[ill, [l6|, the best fit values are = 0.7 and c 2 ^j = — 1 
(compare Figs. 3 and 4). In this case, the LJO pattern 
predicts an overall dynamics (Hubble flow and matter 
fluctuation field) which is for all practical purposes in- 
distinguishable from the traditional A cosmology. Nat- 
urally, such a solution based on a reduction of the dark 
sector, provides not only a possible reinterpretation of 
the ACDM cosmology but also offers a viable alternative 
cosmic scenario to the late time accelerating Universe 
without the need of an exotic dark energy. Indeed one 
of the main advantages of such CCDM cosmology is the 
fact that it contains the same number of free parame- 
ters as the concordance ACDM model, and, therefore, it 
does not require the introduction of any extra fields in 
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its dynamics. 
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